Observation of strong-coupling effects in a diluted magnetic semiconductor Ga^Fe^N 
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A direct observation of the giant Zeeman splitting of the free excitons in Gai-^Fe^N is reported. 
The magnetooptical and magnetization data imply the ferromagnetic sign and a reduced magnitude 
of the effective p-d exchange energy governing the interaction between Fe 3+ ions and holes in GaN, 
-/Vo/3 ( - app ' = +0.5 ± 0.2 eV. This finding corroborates the recent suggestion that the strong p-d 
hybridization specific to nitrides and oxides leads to significant renormalization of the valence band 
exchange splitting. 

PACS numbers: 75.50.Pp, 75.30.Hx, 78.20.Ls, 71.35.Ji 



A strong hybridization between anion p states and 
open d shells of transition metals (TM) is known to ac- 
count for spin-dependent properties of magnetic isola- 
tors, semiconductors, and superconductors, such as an- 
tiferromagnetic superexchange, hole-mediated Zener fer- 
romagnetism, and the contribution of spin fluctuations 
to Cooper pairing. Furthermore, the corresponding ex- 
change splitting of the bands gives rise to giant mag- 
netooptical and magnetotransport phenomena, a finger- 
print of purposeful spintronic materials. Studies of the 
p-d exchange interaction have been particularly reward- 
ing in the case of Mn-based IFVI dilute magnetic semi- 
conductors (DMS) such as (Cd,Mn)Te [H, 3- In these 
systems, Mn is an isoelectronic impurity with a simple 
d 5 configuration, allowing a straightforward, quantitative 
determination of the p-d exchange integral (3 from the so- 
called giant Zeeman effect of the exciton: using the vir- 
tual crystal (VCA) and molecular-field approximations 
(MFA), one calculates a contribution to the spin split- 
tings proportional to the Mn magnetization and to jVq | | , 
where N is the cation density. At the same time, the de- 
termination by photoemission spectroscopy [3| , and the 
computation [J, |5| , of band structure parameters demon- 
strated that the antifcrromagnetic p-d exchange results 
indeed from p-d hybridization. 

According to the above insight, in II- VI oxides and 
III-V nitrides, a small bond length and, thus, strong p- 
d hybridization, should result in particularly large val- 
ues of iVo|/3|, a prediction supported by photoemission 
experiments jfj}. Surprisingly, however, some of the 
present authors found abnormally small exciton split- 
tings in (Zn,Co)0 7], (Zn,Mn)0 8], and (Ga,Mn)N [jj. 
Prompted by this contradiction, one of us suggested [lfj 
that due to the strong p-d coupling, oxides and nitrides 



form an outstanding class of DMS, in which VCA breaks 
down, making the apparent exchange splitting 7Vo/3^ app ^ 
small and of opposite sign. Importantly, the system 
bears some similarity with semiconductor alloys such as 
Ga(As,N), so that experimental and theoretical studies 
of the valence band exchange splitting in the strong cou- 
pling limit may significantly improve our understanding 
of these important alloys. 

To check the above model, we have carried out high- 
resolution studies of magnetization and magnetoreflec- 
tivity in the free exciton region for (Ga,Fe]N epitaxial 
layers, thoroughly characterized previously [111 ]. Since in 
GaN, in contrast to ZnO, the actual ordering of valence 
subbands is settled, the sign of A r o/3^ app ' ) can be unam- 
biguously determined from polarization-resolved magne- 
tooptical spectra. Furthermore, unlike Mn, Fe in GaN 
is an isoelectronic impurity with the simple d 5 configu- 
ration 
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12], allowing a straightforward interpretation 
of the data. Our results lead to a value of iVo/3' app ^ = 
+0.5 ±0.2 eV, which provides an important experimental 
support for the theory [l0| . 

The 0.7 /im thick layers of Ga^Fe^N were grown 
[rl| by metalorganic vapor phase epitaxy (MOVPE), on 
[0001] sapphire substrates with a 1 /im thick, wide-band 
gap (Ga,Al)N buffer layer, which is transparent in the 
free exciton region of Ga^Fe^N. The Fe flow rate was 
adjusted to keep the Fe content well below 0.4%, which is 
the solubility limit of Fe in GaN under our growth condi- 
tions. According to detailed luminescence, electron para- 
magnetic resonance, and magnetic susceptibility studies 
11| . in this range the Fe ions assume mostly the ex- 
pected Fe 3+ charge state corresponding to the d 5 con- 
figuration, for which the spin polarization as a function 
of temperature T and magnetic field B is determined 
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FIG. 1: (color online) (a) Difference between magnetizations 
measured (symbols) at 1.8 and 5 K. Same difference, com- 
puted (solid line) using the Brillouin function for 5 = 5/2 
and treating the Fe concentration x as the only fitting pa- 
rameter, (b) Comparison between the computed magnetiza- 
tion (solid lines, right axis) and the redshift of exciton A in 
a~ polarization (symbols, left axis) for the 0.21% sample at 
three temperatures. The spectra are shown in Fig. . 



by the Brillouin function Bg(T,B) with spin S = 5/2 
and Lande factor gp e 3+ =2.0. This is confirmed by our 
magnetooptical data shown in Fig. [1] which scale with 
B 5 / 2 (T, B). For the reported samples, we determine Fe 
content 0.11±0.02% and 0.21±0.02% by fitting the differ- 
ence in magnetization values measured at 1.8 and 5.0 K 
up to 5 T in a high-field SQUID magnetometer, Fig.QJa). 

Magneto-reflectivity spectra were collected in the Fara- 
day configuration (propagation of light and magnetic 
field along the normal to the sample, which is the c-axis 
of the wurtzite structure), and the light helicity cr^ de- 
fined with respect to the magnetic field direction. Owing 
to the high quality of the layers, either two, or the three, 
free excitons of GaN [ljj], A, B, C, are resolved, and 
their Zeeman shifts are visible in the spectra (Fig. [2]). In 
<7~ polarization, the A — B splitting increases and the 
B — C splitting decreases with the field. Opposite shifts 
are observed in a + . Since these shifts are entirely differ- 
ent from those observed in pure GaN and related 
to the magnetization [Fig. [Tf b)] , we conclude that the 
Fe ions create a "giant Zeeman effect" in (Ga,Fe)N. Re- 
markably, however, the shift observed for A exciton is 
opposite in sign and significantly smaller at given x than 
in other wurtzite DMS with d 5 ions such as (Cd,Mn)Se 

B 

We now describe the methodology that we employed 
to extract the values of the apparent p-d exchange en- 
ergy AT /3( app ) for the valence band and the apparent 
s-d exchange energy iVoa ( - app ) for the conduction band, 
from the reflectivity spectra. It involves two major steps. 
First, we calculate the reflection coefficients as a function 
of the photon energy, using a polariton model which in- 
corporates the ground states of excitons A and B Q, 



but also exciton C, excited states of excitons, and the 
continuum. Second, the field dependence of the exciton 
energies is calculated taking into account the essential 
features of the exciton physics in wide-band gap wurtzite 
DMS UM- 

The starting point of the polariton model is the di- 
electric function e{ui, k) containing polariton poles corre- 
sponding to excitons A and B [7|, [l5j: by solving ana- 
lytically Eq. 3 of Ref. HH, we obtain the refractive index 
and then the reflection coefficient to be compared to ex- 
perimental spectra. However, we replace the background 
dielectric function eg by the residual dielectric function 
€*(u>), which takes into account additional contributions 
which are expected to exhibit no significant polariton ef- 
fect [3]: exciton C, excited states of A, B and C [HI 
which arc optically active (5-states), and transitions to 
the continuum of unbound states [171 ] . Hence 



e*(u>) = Cq + 47ra c 



CJq — lu 2 — iluT. 
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Here eg = 5.2 [lq ]; otQj, luj — Ej/h, and Tj are the polar- 
izability, resonant frequency, and damping rate of each 
exciton A, B, and C, treated as adjustable parameters 
(found to be close to those reported in Ref.lla). Resonant 
energies of the excited states n are huj n j — Ej -\- Rj — 
R*/n 2 , where R* are the effective Rydbergs known from 
our and other [18( studies of excitons in GaN. The cor- 
responding damping parameters are calculated by using 
an empirical formula |l8l.[l9j. r n j = — (r^ — Tj)/n 2 , 
with one common damping rate Too, which is an addi- 
tional fitting parameter. Finally, the contribution £j, u b 
from unbound states is given in Eq. 5 of Ref. [13, with the 
amplitude determined by the exciton polarizabilities aoj 
and the damping parameter Too/2 Calculated 
reflectivity spectra are compared to experimental ones in 
Figs. 2(a,c), and the exciton energies are displayed as a 
function of the magnetic field in Fig. 2(b,d). 

In the second step, the field dependance of the exciton 
energies is calculated for the hamiltonian 



H — E + H v + H c _\ 



H 7 



H 



(app) 
sp—d ' 



(2) 



where Eq is the band-gap energy and H v describes the 
top of the valence band in semiconductors with the 

, ^31. This term includes the 
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wurtzite structure 
trigonal component of both crystal field and biaxial 
strain (described by the parameter Ai) as well as the 
anisotropic spin-orbit interaction (characterized by A2 
and A3 = 5.5 meV [12] for the direction parallel and per- 
pendicular to the c-axis, respectively). The component 
H e -ix describes the electron-hole interaction within the 
exciton BSE El and it involves the effective Rydberg 
and the electron-hole exchange integral 7 = 0.6 meV (23| . 
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FIG. 2: (color online) Reflectivity of Gai-sFe^N in Faraday 
configuration at 1.6 K for x — 0.21% (a) and 0.11% (c) where 
particularly well resolved excitons A, B, and C are visible in 
a~ polarization. (b,d) Field-induced exciton shifts of the exci- 
tons determined from the polariton model (points) compared 
to the expectations of the exciton model (solid lines) . The de- 
termined values of the exciton and exchange parameters are 
shown in Table I. 



Effects linear in the magnetic field are taken into account 
by the standard Zeeman hamiltonian Hz HI, 21 1, param- 
eterized by the effective Lande factor g e = 1.95 for the 
electrons and the relevant Luttinger effective parameter 
k = —0.36 that describes the splitting of all three valence 
subbands [13j. The diamagnetic shift is described by a 
single term quadratic in the magnetic field, H^ nm = dB 2 , 
where d = 1.8 ^eV/T 2 (TJ 



Finally, H^ 9 } is the ex- 



change interaction between Fe ions and carriers in the 
extended states from which the excitons are formed. We 
use the standard form of the s,p-d hamiltonian [U I2I [lil] . 
which is proportional to the scalar product of the car- 
rier spin and the magnetization. However, being aware 
that these splittings can have a different meaning from 
the usual one 1^, 24 1, we use effective quantities a^ app ' 
and /3( app ). The field dependence of the exciton energies, 
calculated with the values of fitting parameters collected 
in Table I, is compared with the experimental results in 
Fig. 2(b,d). 

The detailed analysis of this hamiltonian [9j shows that 



0.11% 
0.21% 



0.44 ± 0.2 
0.46 ±0.1 



0.60 ±0.2 
0.63 ±0.3 



3.4864 19.4 7.5 
3.4846 15.5 7.2 



TABLE I: Fitting parameters describing the exciton energies 
and the giant Zeeman effect shown in Fig. (Hb,d). Units are 
eV, except for Ai and A2, which are given in meV. 



the giant Zeeman shifts of excitons A and B are mainly 
governed by N /3 (app) -AW app) . Moreover, exciton A, 
which in GaN is formed from valence band states with 
parallel spin and orbit (Fg state), shifts to low energies 
in c~ polarization: hence the exchange integral differ- 
ence Af /3 (app )--/Voa (app ) is positive. At the same time, 
exciton B is mixed with exciton C, whose shift is pri- 
marily controlled by iV /3 (app) ±iVoa (app) . Hence, by in- 
corporating exciton C in our description of the reflec- 
tivity spectra, we can evaluate independently the appar- 
ent exchange energy iVo/3 ( - app ' ) for the valence band and 
the apparent exchange energy A r oct;( app ) for the conduc- 
tion band. This determination is quite accurate in the 
case of the 0.11% sample, for which exciton C is spec- 
trally well resolved, but an estimation is still possible 
in the 0.21% sample through its effect on exciton B. 
From the shifts of B and C, we find that the sign of 
jY 0/ g(app)_(_^r oQ ,(app) j s a i g0 p 0S it j ve] an d ft s value is quite 

close to that of iVo/3 (app) -iVoa( app \ see Table I. Hence, 
|iV a( app )| is much smaller than |iV ^ (app) |, and W /3 (app) 
is positive. More precisely, the complete fit yields the 
values of the exchange energies iVo/3 (app) = ±0.5±0.2eV 
and 7V a( app ) = +0.1 ± 0.2 eV. 

We note that our evaluation of A / oa < - app - ) includes 
within the experimental error the values of Nga = 0.25 ± 
0.06 eV found in early studies of Mn-based II- VI DMS 
[2J. We have no reason to question here the applicabil- 
ity of the standard description of the conduction band 
in (Ga,Fe)N - contrary to (Ga,Mn)As and (Ga,Mn)N, 
for which the magnitudes of iVoa:( app - ) were found to be 
reduced under some experimental conditions [24| . 

For the expected d level arrangement, both the ferro- 
magnetic sign and the small magnitude of the apparent 
p-d exchange energy are surprising. Indeed, in GaN, the 
Fe d 5 /d 6 acceptor-like level resides less than 3 eV above 
the top of the valence band [l2j]- This state, e g [, and 
also the higher lying £2 J. level that can hybridize with 
the valence band states, remain unoccupied in intrinsic 
(Ga,Fe)N. At the same time, no donor-like d 5 /d 4 state 
has been found within the GaN gap. This could be ex- 
pected, as in the TM series a particularly large correla- 
tion energy U separates the d 5 and d 6 shells. Hence, the 
occupied Fe ^2! levels reside within the valence band. 
According to the Schrieffer-Wolf theory, in such a case 
the p-d exchange coupling is antiferromagnetic: this was 
confirmed by magnetooptical studies of tellurides and 
selenides containing either Mn or Fe, including stud- 
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ies carried out in our labs, which lead systematically to 
N /3= -1.4 ±0.5 eV @. 

However, it has been recently remarked [l(| that for an 
appropriately strong TM potential, like the one expected 
for oxides and nitrides, the TM ion can bind a hole - a 
trend which was already suggested by strong deviations 
from the VCA in (CdJVIn)S [25j and by the analysis of ab- 
initio calculations [261 ] . A summation of infinite series of 
relevant self-energy diagrams demonstrates that in such 
a situation, the spin splitting of extended states involved 
in the optical transitions remains proportional to magne- 
tization of the localized spins, but the apparent exchange 
energy becomes significantly renormalized (loj . In fact, 
for the expected coupling strength, the theory predicts 
-1 < /3( a PP)//3 < 0, as observed here for (Ga,Fe)N. 

A fruitful comparison can be made with the modifica- 
tion of the conduction band of GaAs induced by a slight 
doping with Nitrogen 27j . When a hydrostatic pressure 
is applied, the Nitrogen isoelectronic centers create local- 
ized states, which are observed in photoluminescence, but 
also extended states (the so-called E + states) to which 
the oscillator strength is transferred so that they are ob- 
served in reflectivity; moreover, these optically active 
states exhibit an anticrossing with the localized states, 
and the strength of the anticrossing increases with the 
Nitrogen density. As a result, the transition to these E + 
states exhibits a shift to high energy when the density of 
low energy states increases. In a DMS with a large value 
of JVo | /3 1 , only the TM impurities with the right spin ori- 
entation (antiparallel to the hole spin) are expected to 
form a localizing center 1^, 25 1: hence the giant Zeeman 
shift in (Ga,Fe)N can be understood as resulting from a 
similar anticrossing but in this case the density of rel- 
evant localized centers is additionally controlled by the 
field-induced orientation of the localized spins. 

In conclusion, giant Zeeman splitting has been ob- 
served by magnetoreflectivity for the A, B, and C ex- 
citons in Gai-^Fe^N. The spectra are well described by 
the exciton model valid for DMS with the wurtzite struc- 
ture. The determined sign and magnitude of the appar- 
ent p-d exchange energy N P^ app ^ = +0.5 ± 0.2 eV con- 
stitutes an important verification of a recent theory [loj ] , 
which describes the effects of the p-d interaction circum- 
venting the virtual-crystal and molecular-field approxi- 
mations that break down in nitrides and oxides. In these 
systems, TM ions bind holes, precluding in this way the 
occurrence of carrier-mediated ferromagnetism in p-type 
materials. However, at sufficiently high hole densities, an 
insulator-to-metal transition is expected. In the metallic 
phase, many-body screening of local potentials annihi- 
lates bound states. Large spin-splitting and robust fer- 
romagnetism are expected in this regime 28, 29j|. 
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